Dynamical Mean Field Theory of Double Perovskite Ferrimagnets 
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The dynamical mean field method is used to analyze the magnetic transition temperature and 
optical conductivity of a model for the ferrimagnetic double perovskites such as SriFeMoO§. The 
calculated transition temperatures and optical conductivities are found to depend sensitively on the 
band structure. For parameters consistent with local spin density approximation band calculations, 
the computed transition temperatures are lower than observed, and in particular decrease dramat- 
ically as band filling is increased, in contradiction to experiment. Band parameters which would 
increase the transition temperature are identified. 

PACS numbers: 



I. INTRODUCTION 

The calculation of non-zero-temperature and dynami- 
cal properties, such as magnetic transition temperatures 
and conductivities, is a long-standing and difficult prob- 
lem in materials theory, but one for which the recent 
theoretical development of the 'dynamical mean field' 
methodcl offers a promising avenue for progress. This 
method, which allows an exact (in principle) treatment 
of local correlations has been recently used in combina- 
tion with 'ab initio' band calculations Jto estimate the 
transition temperatures of Fe and NiU , and in com- 
bination with a tight binding parametrization of band 
theory to elucidate the ptesws of the 'colossal' mag- 
netoresistance manganitesBflCI'El and models of magnetic 
semiconductorslm 

In this paper we apply the method to study ferrimag- 
netic members of the 'double perovskite' family of materi- 
als. These are compounds of chemical formula ABB'Oq, 
with A an alkaline earth such as Sr, Ca or Ba, and 
B, B' two different transition metal ions. In the fer- 
rimagnets of present interest the B site is occupied by 
Fe and B' site by a member of the Ad or 5d transition 
metal series such as Re or Mo. The double perovskite 
family of materials have long been knowno but the fer- 
rimagetic membe.rsJis.ted above have become the subject 
of recent interestE£lE3E30u3 because they seem to be to 
be metallic (except perhaps in the Ca casetJ ferrimag- 
nets with magnetic transition temperatures greater than 
room temperature and apparently highly spin-polarized 
conduction bands, raising the possibility of interesting 
device applications as 'spin valves.'c3, elements in mag- 
netic information storage systemalJ and as soucces of spin 
polarized electrons for spintronic applications^. 

The materials are also of fundamental interest, be- 
cause their physics and materials science are far from 
understood. For example, apparently minor changes 



in chemical composition or processing conditions can 
change the electrical behavior from metallic to insulat- 
ing, or the magnetic transition temperature from 4O0K 
to 200K. More generally, the materials provide examples 
of novel behaviors associated with partial filling of transi- 
tion metal Ad/5d shells, which have been little, studied in 
comparison to the 3d transition metal seriesE3. Further, 
the specific materials we wish to study are believed to be 
half metallic ferrimagnets, but the obvious strong inter- 
action which could give rise to half metallicity (a Hunds 
coupling on the -Fe) exists on only one of the two sites, 
unlike the situation in the somewhat analogous CMR ma- 
terials. Thus behavior arising from the interplay of mag- 
netic order and carrier motion may be different. Finally, 
the origin of the magnetism is not settled. A natural 
guess is that it is due to the strong Hunds coupling on 
the Fe site, but other interactions have been proposed to 
be importantlij 

In this paper we present a model for the low (less than, 
say 3eV) lying electronic states, which are important for 
transport and magnetism. Our model consists of a tight 
binding description of the bands, which we derive from 
previously published first principles calculationstilEj and 
general arguments, and a local interaction (namely a 
Hunds coupling on the Fe site). The key physical as- 
sumption made in our model is that the magnetism is 
driven by the strong Hunds coupling on the Fe site. 
Important technical issues include the two dimensional 
nature of the underlying band structure and the multi- 
orbital nature of the material. We solve the model in 
the dynamical mean field approximation, and from our 
solution determine the magnetic transition temperature 
and optical conductivity, and attempt to determine the 
general materials aspects which control T c . This paper 
supercedes a previous papeiEJ, in which the model Hamil- 
tonian used did not provide an adequate approximation 
to the underlying band structure. 
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II. MATERIAL AND MODEL 
A. Material 

Double perovskite systems form in the ABB'Oq crystal 
structure which generalizes the ABO3 perovskite struc- 
ture familiar from ferroelectrics, high temperature su- 
perconductors and the 'colossal' magnetoresistance rare 
earth manganites by having two different B site ions. In 
the double perovskite materials of interest here, A is an 
alkaline earth such as Sr, Ca or Ba and the B, B' sites 
form a rocksalt structure, i.e. a simple cubic lattice with 
a doubled unit cell and one sublattice occupied by Fe 
and the other by a transition metal from the 4d or bd 
series such as Mo or Re. The crystal fields and atomic 
energetics are such that the formal valences correspond 
to Fe with a half filled, maximally polarized c£-shell while 
the Mo/ Re has,aiie or two d electrons distributed over 
the ti g levelaa'EJO We will focus on electronic states 
arising from the transition metal d-levels. 

B. Hamiltonian 

1. Overview 

The Hamiltonian describing the low lying, electroni- 
cally active degrees of freedom may be written as the 
sum of a 'hopping' part arising from the band structure 
and an interaction part: 

H = Hband + Hi n t (1) 

The relevant portions of the calculated0lil band struc- 
ture involve three bands (degenerate in the ideal double 
perovskite structure) arising from the three transition 
metal t^g levels d xyiyz>xz To a high degree of accuracy 
these three bands do not hybridize with each other and 
the physics is therefore described by a three-fold degen- 
erate tight binding model. The planar character of the 
tig levels implies that the tight binding model has an in- 
teresting two dimensionality, which may be summarized 
as follows. The d xy orbital on a Fe site hybridizes via a 
matrix element t\ with the d xy levels on the four near- 
est neighbor (Mo /Re), sites in the same plane and via a 
much smaller matrix element £3 to the four nearest Fe 
ions also in the same plane. The hopping in the third di- 
rection is negligible, because of the planar character and 
xy orbital symmetry of the d xy wave function. The d xy 
orbital on a Mo/ Re site hybridizes with the four in-plane 
near neighbor Fe sites via the same hopping matrix ele- 
ment t\ and with the four in-plane second nearest neigh- 
bor (Mo/ Re) sites, via another matrix element t 2 , which 
is not particularly small, because of the more spatially 
extended character of the ci-electrons in Ad/5d orbitals. 
Further neighbor hoppings are also found to be impor- 
tant in other *2g-based Ad systems such as Sr2RuOj£Q. 

It is natural to assume that the magnetic character of 
the material derives from the strongly magentic nature 



of the Fe ion and we therefore assume that the dominant 
interaction arises from the strong atomic Hunds coupling 
of the Fe. 



2. Hopping Hamiltonian 

To write the Hamiltonian explicitly we focus the cu- 
bic lattice of B,B' sites in the underlying single per- 
ovskite structure, labelling these sites by i and the op- 
erator creating an electron into orbital a(= xy,yz,xz) 
and spin a by c/ li a . Although we refer to this orbital 
as a 'd-orbital' it in fact represents a hybrid, composed 
mainly of transition metal d and oxygen p orbitals, of the 
correct local symmetry. We introduce a nearest neigh- 
bor (Fe <-> Mo/ Re) hopping t\ and two second neigh- 
bor (same sublattice) hoppings ti and t% representing 
Mo — Mo or Fe — Fe hoppings respectively. As noted 
above we expect that t 2 corresponding to to Mo — Mo 
hopping is relatively large, because of the larger spatial 
extent of the 4d/5d orbitals while £3 is essentially negligi- 
ble. To obtain the conductivity we couple in the electric 
field by using a vector potential and the Peierls phase 
ansatz. This approximation has been shown. to be ac- 
curate in other transition metal oxide contextsE3E3.Thus 
the hopping portion of the Hamiltonian is the sum of 
three identical tight binding models, one for each orbital. 
The Hamiltonians take the general form (note that the 
first sum runs over all lattice sites, the second over the 
B' (non-Fe) and the third over the B (Fe) sites, while 
S a labels the in-plane direction relevant to orbital a and 
we have set the electric charge e and the speed of light c 
equal to unity) 

Hband = - £ ( t ^ eiASaC t,i ,<r C °;i+Sa,<T + H -C.) 

- £ {t 2 ^ A <e+ a ^ +s ^+H.e) (2) 

a,ieB',S' a 

- £ (t 3 , a e* A <c+ aa c l+s ,^ + Hc) 

Hband implies an interesting band structure, which is 
most plainly revealed by writing Hband in momentum 
space in a matrix notation where the upper left entry 
corresponds to Fe and the lower right to Mo, thus if 
A = we have, for the xy orbitals 

Hband,x y [A = 0] = (3) 

( —2ti (cos p x 4- cosp y ) \ 

[ — 2<i (cos + cosp y ) —4t 2 cos p x cosp y J 

where we have set the Fe — Mo distance to unity and 
the momenta are restricted to the reduced Brillouin zone 

\Px\ + \Py \ < T- 
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3. Interaction 

The most important interaction effect constrains the 
occupancy of the B (Fe) site. The formal valence of Fe 
is d 5 and the strong Hunds coupling characteristic of Fe 
(and found in the local spin density approximation to 
band theory) means that in the d 5 configuration all of 
the Fe d-electrons are aligned, leading to a filled, com- 
pletely spin-polarized d-shell. Two charge fluctuation 
processes are possible: Fed 5 «-> Fed 4 or Fed 5 <-> Fed 6 
The strongly stable nature of the filled d-shell implies 
that the d 5 — d 6 process is dominant. To express this 
physics we introduce a strong Hunds coupling on the B 
(Fe) site, expressing the fact that in the ground state the 
Fe is in the d 5 maximal spin configuration, and an energy 
splitting parameter A expressing the differing electroneg- 
ativities of the B and B' sites. Thus we write 

Hint = — J ^ S i ■ C^i^a 'aP c a,i,p + 
a,i£Ba/3 

J2 A'V-,., T r„.,.„ (4) 

For the calculations presented in this paper we will spe- 
cialize to cubic symmetry, so the A are the same for all 
three orbitals, but this restriction may easily be lifted. 
The energy scale relevant for the Fed 5 <-> Fed 6 valence 
fluctuation— is, J — A. Examination of published band 
structureald'Ej indicates that J — A ~ leV while the 
d 5 «-> d A process has a much larger energy of | J + A > 
beV. In a fully spin polarized ground state, the interac- 
tion terms simply become level shifts, A ma j = J + A 
and A m i n = — J + A for the majority and minority spin 
bands. Transitions onto the majority-spin Fe orbital in- 
volve very large energies, so to simplify the calculations 
at T > T c we will take the limit J + A — > oo with J — A 
fixed. We henceforth refer to the quantity J — A as A. 
Because the local spin density approximation may not 
be accurate for strongly interacting systems such as the 
double perovskites, we will consider a range of A here. 

Other authorstj have argued that an additional Hunds- 
type coupling on the B' site is important. Technical lim- 
itations prevent us from treating such an interaction ac- 
curately, so we do not include it here. 

4- Discussion: T=0 Band structure, ferromagnetic case 

Eq|| may be thought of as describing two bands of elec- 
trons: one on the Fe sites, with 'intrinsic' bandwidth set 
by t% and one on the non-Fe sites, with 'intrinsic' band- 
width set by ti. The two bands hybridize via the overlap 
— 2ti (cospz + cosp y ). We see immediately that the hy- 
bridization vanishes along the line cosp K + cosp y = 
which also contains the van Hove points 0, 7r and 7r, at 
which the density of states of the two individual bands 
diverges. Near these points a complicated structure in- 
cluding divergences in the density of states is expected. 




-1.5 -1 -0.5 0.5 1 1.5 
Energy [eV] 

FIG. 1: Total (solid line) and Fe-projected (dashed line) 
densities of states computed using tight binding parameters 
ti = 0.25 eV, t a = 0.15eV, t 3 = 0.03eV. Top panel: T = 
majority spin density of states 5eV Middle 

panel: T = minority spin density of states computed from 
Eq. ^| with A m i n = 0.3eV. Lowest panel: total (both spins) 
density of states at T > T c computed as described in section 
III. 
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The full density of states and the projection of this 
density of states onto the Fe orbitals are shown in the 
upper panels of Fig. 1 for parameters t\ — 0.25 eV, 
£2 = 0.15eV, £3 = 0.03ey (note that most of the ma- 
jority spin Fe density of states is at a low energy out- 
side the range of this plot). Comparison of this (ifp-i 
sity of states to the published band theory resultsy'EHl 
shows that these parameters reproduce the band density 
of states accurately. The main difference is that if the J 
and A are adjusted to correctly reproduce the minority 
spin band, then the upper (non-Fe-portion of the major- 
ity spin band is positioned about 0.5eV too low in energy. 
The extra shift in the majority spin Mo orbitals must be 
attributed to a Hunds coupling on the non-Fe site, not 
included here. 

The two features seen in our calculated density of 
states near 0.5eV arise from states in the vicinity of 
the van Hove points (0, it) and (tt, 0) where the hy- 
bridization vanishes and the B and B' sites have en- 
ergy A + 4^3 ~ OAeV and 42^ ~ 0-6eV respectively, 
whereas the peak at to — arises from the van Hove 
point (it/2, tt/2) of the B' (non-Fe) band, where as noted 
above the hybridization to the Fe vanishes. 

Formal valence arguments indicate that the material 
contains one or two d-electrons beyond the filled shell 
Fe — d 5 Re I Mo — d° configuration. Fig 1 shows that 
for the band theory parameters, these carriers go into 
states with only a small admixture of Fe. The physics 
behind this result is that for this sign of ti the strongly 
hybridized states near p x = p y = are at the bottom of 
the band described by the ti— only term in Hb an d, and 
are pushed further away from the Fe states by usual 
level repulsion, leading to a mainly non-Fe character of 
the lowest states. For n — 1 (Re) only the minority spin 
band is occupied (the majority spins occupy low-lying Fe 
states off of the plotted scale). However, for n = 2 within 
this approximation, chemical potential is \xi « and the 
majority spin band is somewhat occupied, so the mate- 
rial is not a 'half-metal' in this approximation. These 
features will be seen to be of importance for the calcu- 
lated transition temperature and optical conductivity. 

It is interesting to consider a contrasting set of param- 
eters, for which the level repulsion argument works in the 
opposite manner. If £2 has an unphysical (negative) sign 
and A is near then the low lying states of the i2-band 
do not mix with the Fe states, which are pushed down- 
wards by hybridization with the higher- lying Re /Mo lev- 
els, leading to low-lying states of mainly Fe character, as 
shown in Fig 2. 



C. Conductivity 

The current operator J = (5iJ/<5^li!. For electric field 
in the x direction the xy and xz orbitals contribute, thus 
Jx = Jx,xy + Jx.xz with 5 X the lattice vector in the x 
direction and 5' labelling the four 'second neighbor' lat- 
tice vectors (S' = ± (S x ± S y )) so by expanding Eq || in 
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-1.5 -1 -0.5 0.5 1 1.5 
Energy [eV] 

FIG. 2: Total (solid line) and Fe-projected (dashed line) den- 
sities of states computed using alternative tight binding pa- 
rameters U = 0.25 eV, t 2 = -0.15eV, t 3 = 0.03eV. Top 
panel: T = majority spin density of states computed from 
Eq. I using A maj = -2.5eV. Middle panel: T = mi- 
nority spin density of states computed from Eq. |^ with 
A m i„ = 0.3eV. Lowest panel: total (both spins) density of 
states at T > T c computed as described in section III. 
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powers of A we obtain 

Jx,xy(A^) — 

- J2 ( iS xh, a e lA ' s c+ y i a c X y A+5xta - H.c.) - 

(5) 

The expectation value of the term in J proportional to A 
gives the total oscillator strength, S(oo) = f J duj<Ji(uj) 
in the conduction hand contribution to the optical con- 
ductivity (seeErLilS for details). Restoring units we have 
(a is the Fe — Mo/ Re distance, the sum rule is con- 
ventionally defined in terms of the conductivity per unit 
volume and the factor of two comes from the xy and xz 
orbitals, whose contributions to the expectation values 
are identical in cubic symmetry) 



S(oo) 



2e 2 



(6) 



\ieB',S' 



The conductivity is 



S(o )-2 Xjj (n) 



(7) 



with \jj the usual Kubo formula current-current corre- 
lation function evaluated using J X (A = 0) (Eq ^|) and 
again the factor of two represents the orbital degeneracy. 



III. METHOD OF EVALUATION 



G 



xy,band 
BB 



■ UJ 



tt\ (cos(p x ) + cospy) 
u + A*y - At 2 {cos( Px ) cos( Py )) 

(8) 



We measure momenta in units of ir/a where a ks A A is 
the distance from a B to a nearest neighbor B' site. The 
two dimensional Brillouin zone is defined by \p x + p y \ < 

7T. 

The physics is then described by a three-orbital lo- 
cal theory given by the partition function Zi oc = 
J Uc + ce~xp[Si oc ] with an action Si oc which we write in 
the Matsubara frequency representation as 

Sioc - Tj2Tr[c+ a (u) (a^M - JS • a a p) c a p(uj)} (9) 

involving fields c aa and specified by a tensor mean field 
function a which has orbital (ab) and spin [af3) indices 
(the trace is over the spin and orbital indices). The mean 
field function is fixed by the condition that the Green 
function defined from 5; oc , 



G /oc (r) 



loc 



5sl(t) 



(a-S) 



(10) 



is equal to the local Green function defined by integrating 
Eq U over momenta using the self energy defined by Eq 
|gi.e. 



Jpp_ 

(27T) 2 BB 



xy,band 



(11) 



and the integral is over the Brillouin zone defined above. 

Substitution of Eqs ^||l^ into Eq|ll| yields explicit equa- 
tions which arc solved numerically by iteration. 



A. Overview 



B. Calculation of T c 



To evaluate the properties of H we use the dynami- 
cal mean field mcthodEla This method is extensively de- 
scribed and justified elsewhere, and is relevant here be- 
cause the principal interactions are local. In brief the 
central approximation is that the electron self energy, S, 
is momentum independent. In this circumstance the 
physics may be derived from a local theory which may 
be viewed as a quantum impurity model combined with 
a self consistency condition. The application to the dou- 
ble perovskite systems requires some discussion. In these 
systems the unit cell contains two sites, each site con- 
tains three orbitals and there are two choices of spin, so 
the local theory has twelve degrees of freedom. However, 
the problem may be simplified. First, the three orbitals 
(d xy etc) are coupled only via the interaction. Second, 
the interaction exists only on the Fe site, so that we may 
formally integrate out the electrons on the non-Fe (B') 
site, defining a single-orbital model with the inverse F e 
((-B)-site Green function for e.g. the xy orbitals viz 



We calculate the ferromagnetic transition temperature 
by decomposing the mean field function a into non- 
magnetic (do) and magnetic (a\) parts 



afm • a 



(12) 



and linearizing in a\. We take the magnetization direc- 
tion m to be parallel to z and take the limit J — » oo 
so that after integrating out the fermions and redefining 
a — * a + J we obtain (cos(0) is the dot product between 
the direction of the core spin and of the magnetization) 



Tr In [ag(w) + <(w) cos(6»)] 



(13) 



where the Tr is over the frequency index and the orbital 
degree of freedom. 

The Green function of the impurity model becomes 



C? mp (w) 



1 / l-S-a 

2 \ ao — ai cos(0) 



(14) 
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where the expectation value is over the directions of the 
'core spin' S. 

In the paramagnetic phase a\ = 0. Expanding near 
the magnetic transition (assumed second order) yields 



with m =< cos 9 >so that 



S(w) = -ao - ( 2a m - — ) a z 



a i 



(15) 



(16) 



At T > T c m = oi = and substitution of Eqs |T|, |T| 
into Eq [ll] yields 



2a (w) 



where the n order integral /„ is given by 



cPp 



(17) 



(18) 



This equation is solved numerically by iteration for a 
sufficiently dense set of frequency points (typically fre- 
quency spacing 0.04ti). Once a solution for ao is obtained 
we may linearize Eq O in the magnetic part of the self 
energy and local Green function, obtaining 



m 
2a 



«i 
6a 2 Q 



i(u,a (co)) (2a Q m 



ai 
3 



where 



m =< cos >= 



(19) 



(20) 



for n = 2 (Mo case) the ground state is not ferromag- 
netic. These calculated values are in sharp disagreement 
with the experimental values T c ^ A00K for both n = 1 
and n = 2. The relatively small values of T c found for 
n = 1 in this calculation may be understood from the 
density of states, which shows that the low-lying states 
lie mainly on the non — Fe sites, which are far displaced 
in energy from the magnetic site and therefore do not 
hybridize strongly with it, so the effective carrier-spin in- 
teraction is not strong. That the n = 2 is non-magnetic 
may be understood by combining the results ofl with 
the observation that the band structure is effectively two 
dimensional. In the extreme weak coupling limit, the 
nature of the magnetic ground state is determined by 
the wave vector at which the susceptibility is maximal. 
For the two dimensional band structures considered here 
this maximum is not at q = 0. Refu showed that in the 
DMFT approximation, increasing the carrier-spin cou- 
pling increased the range in which ferromagnetism ex- 
isted, but that as band filling is increased, a transition 
to an antiferromagnetic state generically occurs, and gets 
pushed to the half-filled band only for J of the order of 
the bandwidth. These effects are more pronounced for 
the two dimensional band structure we consider. 

The relative weakness of the virtual Re /Mo <-» Fe 
transitions is reflected in the temperature dependence of 
the many-body density of states, shown for the 'band' 
parameters in the lower panel of Fig. 1. Comparing 
these we see that disordering the Fe spins leads to a 
slight narrowing of the bands, but the larger (3|Q%) band 
narrowing effects found in CMR manganitesu are not 
observed for these parameters. 



Solving for a\ and then using this to obtain an expression 
for m yields a self consistent equation for T c which in the 
limit of cubic symmetry becomes 



^ 1 + 2alh ^ 



1 + 2a^I 2 



(21) 



It turns out that the transition temperatures are low 
compared to the other scales of the model so that that 
one may recast this equation via analytical continuation 
to the real axis as (/i is the chemical potential correspond- 
ing to the desired carrier density) 



T = l" — Im 

, 7T 
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1 + 24h 



(22) 



IV. RESULTS AND DISCUSSION 



n=l 



n=2 





FIG. 3: Left panel: Contour plots of calculated transtion tem- 
peratures for range of model parameters and conduction band 
density n = 1. Right panel: Contour plots of calculated tran- 
sition temperature for n — 2. In each figure, contours are 
spaced approximately 100 Kelvin apart and the white areas 
correspond to transition temperatures in excess of 400K. 



We have calculated the magnetic transiton tempera- 
ture from Eq. ^2|, finding that for the parameters used 
to construct Fig. 1 (and which are the ones following 
from band theory) n = 1 (Re case) T c w 110K and that 



To understand the behavior of the model in more detail 
we have evaluated the predicted ferromagnetic transtion 
temperatures for wide range of model parameters. The 
results are summarized in the two panels of Fig. 3 which 
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show via contour plots the values of T c predicted by the 
method. The contours are spaced approximately 100A" 
apart, and the black regions indicate the areas in which 
the calculated T c vanishes. It is seen that in order to ob- 
tain a reasonably high transition temperature, especially 
for the n = 2 band filling, one must choose the parameter 
ti to have the opposite sign from the physical one. The 
reason for this behavior is reveal by Fig. 2, which shows 
the density of states for parameters which maximize the 
n = 2 T c . The low- lying states for this case are seen to 
be of mainly Fe character, because the level repulsion 
argument which pushed down the non-Fe states for the 
LSDA parameters is not operative here. 



30000 



c 




o 1 1 

0.25 0.5 0.75 1 

co(eV) 

FIG. 4: T = (dashed) and T > T c conductivities for n = 1 
using best-fit band parameters ti = 0.25eV, t% = 0.15eF 
U = 0.03eF and A = 0.3eF used in Fig. 1. The T = 
conductivity was computed using an artificial broadening of 
Q.leV applied to the B site. 

We have also calculated the optical conductivity for 
various model parameters. Results obtained using the 
LSDA parameters are shown in Figs. 4 (i?e-case, n = 1) 
and fig. 5 (Afo-case, n = 2). One would in principle ex- 
pect two classes of transitions: a 'Drude' peak centered 
at lo — involving motion of electrons near the Fermi 
surface and an interband transition involving moving an 
electron from a Re/ Mo to a Fe. Our calculations indi- 
cate that for the tight binding parameters corresponding 
to the LSDA calculation, the interband feature is very 
weak, indeed not visible in the Figure again demonstrat- 
ing the weakness of the Fe — Mo coupling for these pa- 
rameters. We observe that the 'Drude' part has a dis- 
tinctly non- Drude form, which arises because in our cal- 
culation the scattering processes couple to the 'B' (Fe) 
site only; although the regions of momentum space where 
the hybridization vanishes are of measure zero, they do 
lead to a frequency dependence of the scattering rate 
which explains the peculiar form. We also note that the 
main cause of the changes in conductivity and oscillator 
strength between T = and T = T c is the change in 
band filling, which leads to a change in optical matrix 



30000 




0.5 1 1.5 

(0(eV) 

FIG. 5: Upper panel: T = (dashed) and T > T c (solid) 
conductivities for n — 2. using best-fit band parameters 
ti = 0.25eV, ta = 0.15eF t 3 = 0.03eF and A = 0.3eF used 
in Fig. 1. The T = conductivity was computed using an 
artificial broadening of O.leV applied to the B site. The small 
contribution to the T = conductivity arising from the mi- 
nority spin band is shown as the dash-dot line. 



element. In the paramagnetic state one has three bands, 
each with a two-fold spin degeneracy, corresponding to a 
filling of n/6 (n=l,2 is the particle density), so the rel- 
evant states are quite close to the bottom of the band 
where the optical matrix element is small. In the ferri- 
magnetic state for n = 1 one loses the spin degeneracy, so 
one has three bands each filled to a higher level, so with 
a correspondingly higher fermi velocity and optical ma- 
trix element, whereas for the n — 2 case the temperature 
induced shift corresponds to a change from 1/3 to 2/3 
filled band, with much smaller change in optical matrix 
elements. The temperature dependent change in the os- 
cillator strength is therefore much less. As noted above, 
in the ferrimagentic case for n = 2 one has a small filling 
of the majority spin band, leading to a small additional 
contribution to er, shown as the dot-dashed line in Fig. 
5. 

The conductivity corresponding to the parameters 
which maximize T c (Fig 2) is shown in Fig. 6. We see 
that the different electronic structure leads to a differ- 
ent optical conductivity: the Drude absorbtion is weaker, 
and a peak corresponding to excitation of carriers from 
Fe to Re I Mo is evident. 

The conductivity of SLcaFeMoOe has been measured 
by Jung and co-workersE-l These authors found a con- 
ductivity which was of roughly the Drude form, (albeit 
with a rather larger scattering rate than we have used) 
but additionally has a weak kink at a frequency of ap- 
proximately 0.6eV. It is interesting to speculate that 
this kink is a signature of the 'interband' feature which 
we found only for the 'antiphysical parameters. A more 
detailed experimental investigation of the band structure 
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FIG. 6: T = (dashed) and T > T c (solid) optical con- 
ductivities for n = 1 and band parameters ti = 0.25eV, 
i2 = — O.lSel 7 and A = corresponding to the density of 
states shown in Fig 2. 

may be warranted. 

V. CONCLUSIONS 

We have used the dynamical mean field method to de- 
termine the ferromagnetic transition temperature, den- 
sity of states and optical conductivity of a model repre- 
senting key physics (two dimensionality of band structure 
and strong on-site interaction on Fe site) of the double 
perovskite ferrimagnets Sr2Fe(Mo/ Kc)Oq. Our method 
can easily be generalized to include the effects of mis-site 



disorder, or lattice distortions which split the ti g levels. 
However, such generalization is not immediately war- 
ranted because the calculated transition temperatures 
are, at least for the parameters following from band struc- 
ture calculations, in qualitative disagreement with exper- 
imental data-in particular, the calculation predicts that 
Sr<zFeM.oO% is not ferrimagnetic, whereas experiment in- 
dicates that it is with a T c is excess of A00K , and under- 
predicts the T c of Sr^Fe Re Oq by a factor of almost 4. 

The essential reason for this was found to be that the 
band theory parameters imply that the mobile carriers 
reside mainly on the non — Fe sites, and hybridize weakly 
with these sites. Optical conductivity signatures of the 
weak hybridization were demonstrated. The calculation 
indicates that transition temperatures would be substan- 
tially raised if parameters are used for which the added 
carriers are largely on the Fe sites. An alternative pos- 
sibility is that an interaction omitted from the model 
is crucially important; in particular that the magnetism 
should not be regarded as arising from correlations on 
the Fe site, but should isntead be thought of more as a 
Stoner instability of the band arising from the_Re/Mo 
states. The additional interaction proposed inE£l would 
tend to produce this physics and an important next step 
would be to extend the methods developed here to the 
treatment of this case. 
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